2. The Lagrangian Formalism 



When I was in high school, my physics teacher called me down one day after 
class and said, "You look bored, I want to tell you something interesting" . 
Then he told me something I have always found fascinating. Every time 
the subject comes up I work on it. 

Richard Feynman 

Feynman's teacher told him about the "Principle of Least Action" , one of the most 
profound results in physics. 

2.1 The Principle of Least Action 

Firstly, let's get our notation right. Part of the power of the Lagrangian formulation 
over the Newtonian approach is that it does away with vectors in favour of more general 
coordinates. We start by doing this trivially. Let's rewrite the positions of N particles 
with coordinates as where A — 1,. . . 3N. Then Newton's equations read 

dV 

P^ = -g^ (2.1) 

where = rriAX^. The number of degrees of freedom of the system is said to be 3iV. 
These paramctcrise a 3A^-dimensional space known as the configuration space C. Each 
point in C specifics a configuration of the system (i.e. the positions of all N particles). 
Time evolution gives rise to a curve in C. 



Figure 2: The path of particles in real space (on the left) and in configuration space (on the 
right) . 



The Lagrangian 

Define the Lagrangian to be a function of the positions and the velocities x^ of all 
the particles, given by 

L(a;^, i^) = T{x^) - V{x^) (2.2) 
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where T = ^YlA^^i^^)^ kinetic energy, and V{x^) is the potential energy. 

Note the minus sign between T and V\ To describe the principle of least action, we 
consider all smooth paths x^(t) in C with fixed end points so that 

X^iU) — ^^kiitial ^^(^/) — ^final 



(2.3) 



Of all these possible paths, only one is the true path 
taken by the system. Which one? To each path, let us 
assign a number called the action S defined as 



L{x^{t),x^{t)) dt 



(2.4) 



The action is a functional (i.e. a function of the path which 
is itself a function) . The principle of least action is the fol- 
lowing result: 




Figure 3: 



Theorem (Principle of Leeist Action): The actual path taken by the system is an 
extremum of S. 



Proof: Consider varying a given path slightly, so 

x^{t) ^ x^{t) + Sx^{t) 



(2.5) 



where we fix the end points of the path by demanding Sx^{ti) = Sx^{tf) = 0. Then 
the change in the action is 



ss 



Ldt 



tf 



6L dt 



( ^ 



Sx^ + ^5x^ 1 dt 
ox^ 



At this point we integrate the second term by parts to get 



dL d f dL\ 



dx^ dt \dx^ J 



Sx"^ dt + 



dL 



5x' 



(2.6) 



(2.7) 



But the final term vanishes since we have fixed the end points of the path so Sx'^iti) = 
Sx^{tf) = 0. The requirement that the action is an extremum says that 6S = for all 
changes in the path 5x'^{t). We see that this holds if and only if 



dL d f dL 



dx^ dt V <9i^ 



= foreach A = 1,...3A^ 



(2.8) 
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These are known as Lagrange's equations (or sometimes as the Euler-Lagrange equa- 
tions). To finish the proof, we need only show that Lagrange's equations are equivalent 
to Newton's. From the definition of the Lagrangian (2.2), we have 9L/9a;^ = —dV/dx^, 
while dL/dx^ — pa- It's then easy to see that equations (2.8) are indeed equivalent to 



Some remarks on this important result: 

• This is an example of a variational principle which you already met in the epony- 
mous "variational principles" course. 

• The principle of least action is a slight misnomer. The proof only requires that 
5S — 0, and does not specify whether it is a maxima or minima of 5". Since 
L — T — V, we can always increase S by taking a very fast, wiggly path with 
T ^ 0, so the true path is never a maximum. However, it may be either a 
minimum or a saddle point. So "Principle of stationary action" would be a more 
accurate, but less catchy, name. It is sometimes called "Hamilton's principle" . 

• All the fundamental laws of physics can be written in terms of an action principle. 
This includes electromagnetism, general relativity, the standard model of particle 
physics, and attempts to go beyond the known laws of physics such as string 
theory. For example, (nearly) everything we know about the universe is captured 
in the Lagrangian 



where the terms carry the names of Einstein, Maxwell (or Yang and Mills) and 
Dirac respectively, and describe gravity, the forces of nature (electromagnetism 
and the nuclear forces) and the dynamics of particles hke electrons and quarks. If 
you want to understand what the terms in this equation really mean, then hang 
around for Part III next year! 

• There is a beautiful generalisation of the action principle to quantum mechan- 
ics due to Feynman in which the particle takes all paths with some probability 
determined by S. We will describe this in Section 4.8. 

• Back to classical mechanics, there are two very important reasons for working with 
Lagrange's equations rather than Newton's. The first is that Lagrange's equations 
hold in any coordinate system, while Newton's are restricted to an inertial frame. 
The second is the ease with which we can deal with constraints in the Lagrangian 
system. We'll look at these two aspects in the next two subsections. 



(2.1). 



□ 



(2.9) 
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2.2 Changing Coordinate Systems 

We shall now show that Lagrange's equations hold in any coordinate system. In fact, 
this follows immediately from the action principle, which is a statement about paths 
and not about coordinates. But here we shall be a little more pedestrian in order to 
explain exactly what we mean by changing coordinates, and why it's useful. Let 

Qa^ Qaixi,. . . ,X3N:t) (2.10) 

where we've included the possibility of using a coordinate system which changes with 
time t. Then, by the chain rule, we can write 

dQa dqa .A , dqa , . 

Qn = = — tX + — — (2.11) 

dt dx^ dt ^ ' 

In this equation, and for the rest of this course, we're using the "summation convention" 
in which repeated indices are summed over. Note also that we won't be too careful 
about whether indices are up or down - it won't matter for the purposes of this course. 
To be a good coordinate system, we should be able to invert the relationship so that 
_ 2;^(qr^^ t) which we can do as long as we have det{dx'^ / dqa) ^ 0. Then we have, 

. A dx^ . dx^ . ^ ^, 

Now we can examine L{x^^x^) when we substitute in a;^(g(j,i). Using (2.12) we have 

dl^_dl^dx^ dL / d'^x^ . 

dqa dx^ dqa dx^ \dqadq},^^ dtdqa) 

while 

dL dL dx^ 



dqa dx^ dqa 



(2.14) 



We now use the fact that we can "cancel the dots" and dx^/dqa = dx^/dqa which we 
can prove by substituting the expression for x^ into the LHS. Taking the time derivative 
of (2.14) gives us 



d f dL\ d f dL\ dx^ dL f d'^x^ . d^x^ 



dt \dqaj dt \dx^ J dqa dx^ \dqadqb dqadt 
So combining (2.13) with (2.15) we find 



d_ (dL \ _ dL_ 
dt \dqa) dqa 



dL 



d ( dL\ 

dt \dx^ J dx^ 



dx^ 
dqa 



(2.16) 
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Equation (2.16) is our final result. We see that if Lagrange's equation is solved in the 
x"^ coordinate system (so that [. . .] on the RHS vanishes) then it is also solved in the 
qa coordinate system. (Conversely, if it is satisfied in the coordinate system, so the 
LHS vanishes, then it is also satisfied in the x-^ coordinate system as long as our choice 
of coordinates is invertible: i.e (\et{dx-^ / dqa) ^ 0). 

So the form of Lagrange's equations holds in any coordinate system. This is in 
contrast to Newton's equations which are only valid in an inertial frame. Let's illustrate 
the power of this fact with a couple of simple examples 

2.2.1 Example: Rotating Coordinate Systems 

Consider a free particle with Lagrangian given by 

L = \mv^ (2.17) 

with r = (x, y, z). Now measure the motion of the particle with respect to a coordinate 
system which is rotating with angular velocity a; = (0,0, tu) about the z axis. If 
r' = (x', 7/', z'^ are the coordinates in the rotating system, we have the relationship 

x' — X cos cut + y sin cut 
y' — y cos cut — X sin cut 

z! ^ z (2.18) 

Then we can substitute these expressions into the Lagrangian to find L in terms of the 
rotating coordinates, 

L = \m{[x' - cuy'f + {y' + cux'f + z^\ = |m(r' + w x r')^ (2.19) 

In this rotating frame, we can use Lagrange's equations to derive the equations of 
motion. Taking derivatives, we have 

^ = mix' X UJ - a; X (w X r')) 
or' 

so Lagrange's equation reads 

Jt \di') ~d?^ "^^"^ + u; X (a; X r') + 2u; X r') = (2.21) 

The second and third terms in this expression are the centrifugal and coriolis forces 
respectively. These are examples of the "fictitious forces" that you were warned about in 
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the first year. They're called fictitious because they're a consequence of the reference 
frame, rather than any interaction. But don't underestimate their importance just 
because they're "fictitious" ! According to Einstein's theory of general relativity, the 
force of gravity is on the same footing as these fictitious forces. 

The centrifugal force Fccnt = —muj x (cj x r') points outwards in the plane perpen- 
dicular to (jj with magnitude ma;^|r^| = m|vj_p/ |r'_i_| where j_ denotes the projection 
perpendicular to uj. 




Figure 4: In the northern hemisphere, a particle is deflected in a clockwise direction; in the 
southern hemisphere in an anti-clockwise direction. 

The coriolis force Fcor = — 2ma; x r' is responsible for the large scale circulation of 
oceans and the atmosphere. For a particle travelling on the surface of the rotating 
earth, the direction of the coriolis force is drawn in figure 4. We see that a particle 
thrown in the northern hemisphere will be seen to rotate in a clockwise direction; a 
particle thrown in the southern hemisphere rotates in an anti-clockwise direction. For 
a particle moving along the equator, the coriolis force points directly upwards, so has 
no effect on the particle. 

Questions 

• The coriolis force is responsible for the formation of hurricanes. These rotate 
in different directions in the northern and southern hemisphere, and never form 
within 500 miles of the equator where the coriolis force is irrelevant. But hur- 
ricanes rotate anti- clockwise in the northern hemisphere. This is the opposite 
direction from what we deduced above for a projected particle! What did we 
miss? 

• The Rhine fiows roughly north and straight. Which bank has its river defences 
repaired most frequently? 
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• Estimate the magnitude of the coriohs force. Do you think that it really affects 
the motion of water going down a plughole? What about the direction in which 
a CD spins? 

• Stand on top of a tower at the equator and drop a ball. As the ball falls, the 
earth turns underneath from west to east. Does the ball land 

1 . At the base of the tower? 

2. To the east? 

3. To the west? 

2.2.2 Example: Hyperbolic Coordinates 

A particle moves in the {x, y) plane with a force 
directed towards the origin O with magnitude propor- 
tional to the distance from O. How does it move? In 
Cartesian coordinates, this problem is easy. We have 
the Lagrangian 

L = lm{x'' + f) - ^k{x^ + y^) (2.22) 

Let's set m = k = 1 for simplicity. The equation of 
motion for this system is simply 

X = —X and ij = —y (2.23) 

Now suppose we want to know the motion of the system in hyperbolic coordinates 
defined as 

2xy = II , -y^ = X (2.24) 

The coordinates /i and A are curvilinear and orthogonal (i.e. two hyperbolics intersect 
at 90°). We could try solving this problem by substituting the change of coordinates 
directly into the equations of motion. It's a mess. (Try if you don't believe me!). 
A much simpler way is to derive expressions for x, y, x and y in terms of the new 
coordinates and substitute into the Lagrangian to find, 

i^=l4=^-iAA^ (2-25) 




Figure 5: Hyperbolic coordi- 
nates. 
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Prom which we can easily derive the equation of motion for A 

d fdL\ dL _ d ( A \ .2. A i A 

= (2.26) 

Which is also a mess! But it's a mess that was much simpler to derive. Moreover, we 
don't need to do any more work to get the second equation for fi: the symmetry of the 
Lagrangian means that it must be the same as (2.26) with A interchanged. 

2.3 Constraints and Generalised Coordinates 

Now we turn to the second advantage of the Lagrangian formulation. In writing — 
—ViV, we implicitly assume that each particle can happily roam anywhere in space 
R^. What if there are constraints? In Newtonian mechanics, we introduce "constraint 
forces" . These are things hke the tension of ropes, and normal forces apphed by surfaces. 
In the Lagrangian formulation, we don't have to worry about such things. In this 
section, we'll show why. 

An Example: The Pendulum 

The simple pendulum has a single dynamical degree of freedom 
9, the angle the pendulum makes with the vertical. The position of 
the mass m in the plane is described by two cartesian coordinates x 
and y subject to a constraint x'^+y^ — P. We can parameterise this 
ets X — IsinO and y — I cos 6. Employing the Newtonian method 
to solve this system, we introduce the tension T as shown in the 
diagram and resolve the force vectors to find, 

mx — —Tx/l , my = mg — Ty/l (2.27) 

To determine the motion of the system, we impose the constraints 
at the level of the equation of motion, and then easily find 

e^-{g/l)sme , T ^mW^ + mg cose (2.28) 

While this example was pretty straightforward to solve using Newtonian methods, 
things get rapidly harder when we consider more complicated constraints (and we'll 
see plenty presently). Moreover, you may have noticed that half of the work of the 
calculation went into computing the tension T. On occasion we'll be interested in this. 
(For example, we might want to know how fast we can spin the pendulum before it 




V mg 

Figure 6: 
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breaks). But often we won't care about these constraint forces, but will only want to 
know the motion of the pendulum itself. In this case it seems like a waste of effort to 
go through the motions of computing T. We'll now see how we can avoid this extra 
work in the Lagrangian formulation. Firstly, let's define what we mean by constraints 
more rigorously. 

2.3.1 Holonomic Constraints 

Holonomic Constraints are relationships between the coordinates of the form 

fa{xA,t)^0 a^l,...,3N-n (2.29) 

In general the constraints can be time dependent and our notation above allows for 
this. Holonomic constraints can be solved in terms of n generalised coordinates qi, 
i — 1, . . .n. So 

XA^XA{qi,---,qn) (2.30) 

The system is said to have n degrees of freedom. For the pendulum example above, 
the system has a single degree of freedom, q = 9. 

Now let's see how the Lagrangian formulation deals with constraints of this form. 
We introduce 3N — n new variables Aq, called Lagrange multipliers and define a new 
Lagrangian 

L' = L{x^,x^)^KUx'',t) (2.31) 

We treat Aq hke new coordinates. Since L' doesn't depend on Aq, Lagrange's equations 
for Aq, are 



Fit' 

^ = A(x^i) = (2.32) 

which gives us back the constraints. Meanwhile, the equations for are 

d { dL\ dL _ dfa , . 



dt \dx^ J dx^ " dx^ 

The LHS is the equation of motion for the unconstrained system. The RHS is the 
manifestation of the constraint forces in the system. We can now solve these equations 
as we did in the Newtonian formulation. 
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The Pendulum Example Again 

The Lagrangian for the pendulum is given by that for a free particle moving in the 
plane, augmented by the Lagrange multiplier term for the constraints. It is 

L' = \m{x^ + f ) + mgy + \\{x^ + - t^) (2.34) 

Prom which we can calculate the two equations of motion for x and y, 

mx = Xx and ij = mg + Xy (2.35) 

while the equation of motion for A reproduces the constraint x'^+y'^—l'^ — 0. Comparing 
with the Newtonian approach (2.27), we again see that the Lagrange multiplier A is 
proportional to the tension: A = —T/l. 

So we see that we can easily incorporate constraint forces into the Lagrangian setup 
using Lagrange multipliers. But the big news is that we don't have to! Often we don't 
care about the tension T or other constraint forces, but only want to know what the 
generalised coordinates Qi are doing. In this case we have the following useful theorem 

Theorem: For constrained systems, we may derive the equations of motion directly 
in generalised coordinates qi 

L[qi, Qi, t] = L[x^(g,, t),x^(qi, Qi, t)] (2.36) 
Proof: Let's work with L' — L + Xafa and change coordinates to 




(2.37) 



n 



We know that Lagrange's equations take the same form in these new coordinates. In 
particular, we may look at the equations for q^, 

dt \dqij dqi " % 

But, by definition, d fa/ dqi — 0. So we are left with Lagrange's equations purely in 
terms of qi, with no sign of the constraint forces. If we are only interested in the 
dynamics of the generalised coordinates qi, we may ignore the Lagrange multipliers 
and work entirely with the unconstrained Lagrangian L{qi, qi, t) defined in (2.36) where 
we just substitute in xa — XA{qi, t). □ 



- 19 - 



The Pendulum Example for the Last Time 

Let's see how this works in the simple example of the pendulum. We can parameterise 
the constraints in terms of the generalised coordinate 9 so that x = I sin 6 and y = 
I cos 9. We now substitute this directly into the Lagrangian for a particle moving in 
the plane under the effect of gravity, to get 



From which we may derive Lagrange's equations using the coordinate 6 directly 



which indeed reproduces the equation of motion for the pendulum (2.28). Note that, 
as promised, we haven't calculated the tension T using this method. This has the 
advantage that we've needed to do less work. If we need to figure out the tension, we 
have to go back to the more laborious Lagrange multiplier method. 

2.3.2 Non-Holonomic Constraints 

For completeness, let's quickly review a couple of non-holonomic constraints. There's 
no general theory to solve systems of this type, although it turns out that both of the 
examples we describe here can be solved with relative ease using different methods. We 
won't discuss non-holonomic constraints for the rest of this course, and include a brief 
description here simply to inform you of the sort of stuff we won't see! 

Inequalities 

Consider a particle moving under gravity on the outside of a sphere of radius R. It is 
constrained to satisfy z^ + y'^ + z'^ > E?. This type of constraint, involving an inequality, 
is non-holonomic. When the particle lies close to the top of the sphere, we know that 
it will remain in contact with the surface and we can treat the constraint effectively as 
holonomic. But at some point the particle will fall off. To determine when this happens 
requires different methods from those above (although it is not particularly difficult). 

Velocity Dependent Constraints 

Constraints of the form (/(x"^, x^, t) = which cannot be integrated to give f{x^, t) = 
are non-holonomic. For example, consider a coin of radius R rolling down a slope as 
shown in figure 7. The coordinates {x, y) fix the coin's position on the slope. But the 
coin has other degrees of freedom as well: the angle 9 it makes with the path of steepest 



L = \m{x^ y^) mgy 
— \ml^9'^ -\- mgl cos 9 



(2.39) 




(2.40) 
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Figure 7: The coin rolling down a slope leads to velocity dependant, non-holonomic con- 
straints. 



descent, and the angle that a marked point on the rim of the coin makes with the 
vertical. If the coin rolls without slipping, then there are constraints on the evolution 
of these coordinates. We must have that the velocity of the rim is Vrim = R4>- So, in 
terms of our four coordinates, we have the constraint 

x = R^sm9 , y = R^ cos 9 (2.41) 

But these cannot be integrated to give constraints of the form f{x, y, 9, 0) = 0. They 
are non-holonomic. 

2.3.3 Summary 

Let's review what we've learnt so far. A system is described by n generalised coordinates 
Qi which define a point in an n-dimensional configuration space C. Time evolution is a 
curve in C governed by the Lagrangian 



such that the Qi obey 



L{qi,qi,t) (2.42) 



These are n coupled 2°^^ order (usually) non-linear differential equations. Before we 
move on, let's take this opportunity to give an important definition. The quantity 

P^-^ (2.44) 
dqi 

is called the generalised momentum conjugate to qi. (It only coincides with the real 
momentum in Cartesian coordinates). We can now rewrite Lagrange's equations (2.43) 
as Pi — dL/dqi- The generalised momenta will play an important role in Section 4. 



- 21 - 



Note: The Lagrangian L is not unique. We may make the transformation 



L' = 



aL 



for a e R 



or 



L' 




(2.45) 



for any function / and the equations of motion remain unchanged. To see that the last 
statement is true, we could either plug L' into Lagrange's equations or, alternatively, 
recall that Lagrange's equations can be derived from an action principle and the ac- 
tion (which is the time integral of the Lagrangian) changes only by a constant under 
the transformation. (As an aside: A system no longer remains invariant under these 
transformations in quantum mechanics. The number a is related to Planck's constant, 
while transformations of the second type lead to rather subtle and interesting effects 
related to the mathematics of topology) . 

2.3.4 Joseph-Louis Lagrange (1736-1813) 

Lagrange^ started off life studying law but changed his mind and turned to mathematics 
after reading a book on optics by Halley (of comet fame). Despite being mostly self- 
taught, by the age of 19 he was a professor in his home town of Turin. 

He stayed in Italy, somewhat secluded, for the next 11 years although he commu- 
nicated often with Euler and, in 1766, moved to Berhn to take up Euler's recently 
vacated position. It was there he did his famous work on mechanics and the calculus of 
variations that we've seen above. In 1787 he moved once again, now to Paris. He was 
just in time for the French revolution and only survived a law ordering the arrest of 
all foreigners after the intervention of the chemist Lavoisier who was a rather powerful 
political figure. (One year later, Lavoisier lost his power, followed quickly by his head.) 

Lagrange published his collected works on mechanics in 1788 in a book called "Mechanique 
Analytique" . He considered the work to be pure mathematics and boasts in the intro- 
duction that it contains no figures, thereby putting the anal in analytique. 

Since I started with a quote about Newton's teaching, I'll include here a comment 
on Lagrange's lectures by one of his more famous students: 

"His voice is very feeble, at least in that he does not become heated; he 
has a very pronounced Italian accent and pronounces the s like z ... The 
students, of whom the majority are incapable of appreciating him, give him 
little welcome, but the professors make amends for it." 



^You can read all about the lives of mathematicians at http://www-gap.dcs.st-and.ac.uk/ his- 
tory/Bioglndex.html 



Fourier analysis of Lagrange 
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2.4 Noether's Theorem and Symmetries 

In this subsection wc shall discuss the appearance of conservation laws in the Lagrangian 
formulation and, in particular, a beautiful and important theorem due to Noether 
relating conserved quantities to symmetries. 

Let's start with a definition. A function F{qi,qi,t) of the coordinates, their time 
derivatives and (possibly) time t is called a constant of motion (or a conserved quantity) 
if the total time derivative vanishes 

dF ^ fdF OF \ OF , ^ , 

whenever qi{t) satisfy Lagrange's equations. This means that F remains constant along 
the path followed by the system. Here's a couple of examples: 

Claim: If L does not depend explicitly on time t (i.e. dL/dt — 0) then 

is constant. When H is written as a function of qi and pi, it is known as the Hamiltonian. 
It is usually identified with the total energy of the system. 

Proof 

dH ^ f dL d fdL\ dL dL \ , , 

^ = E l^^^a^ + 'iut ) - - ^'-''^ 

which vanishes whenever Lagrange's equations (2.43) hold. □ 

Claim: Suppose dL/dqj — for some qj. Then qj is said to be ignorable (or cyclic). 
We have the conserved quantity 

ft = ^ (2.49) 

Proof: 

dpj d f dL\ dL 



dt dt \dqj J dqj ^ (2.50) 
where we have used Lagrange's equations (2.43) in the second equality. □ 
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2.4.1 Noether's Theorem 

Consider a one-parameter family of maps 

qi{t) Qi{s,t) 



seR 



(2.51) 



such that Qi{0, t) = qi(t). Then this transformation is said to be a continuous symmetry 
of the Lagrangian L if 



d_ 
ds 



L{Qi{s,t),Qi{s,t),t) = 



(2.52) 



Noether's theorem states that for each such symmetry there exists a conserved quantity. 
Proof of Noether's Theorem: 



dL dL dQi ^ dL dQi 



ds dQi ds dQi ds 



(2.53) 



so we have 



dL 
ds 



_ dL dQ, 

s=0 ^« 



dL dQ, 



s=0 



d fdL\ dQi 



dt \ dqi ) ds 

_ d_ (dL dQi 



dqi ds 
+ 



s=0 

dL dQi 



dqi ds 



(By Lagrange) 



dt \dqi ds 



s=0 



(2.54) 



and the quantity Yl-{dL/ dqi){dQi/ ds) , evaluated at s = 0, is constant for all time. □ 



Example: Homogeneity of Space 

Consider the closed system of N particles discussed in Section 1 with Lagrangian 



(2.55) 



This Lagrangian has the symmetry of translation: r^ — >■ rj + sn for any vector n and 
for any real number s. This means that 



L{ri, Ti, t) = L{ri + sn, r^, t) 



(2.56) 



This is the statement that space is homogeneous and a translation of the system by 
sn does nothing to the equations of motion. These translations are elements of the 
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Galilean group that wc met in section 1.2. From Noether's theorem, we can compute 
the conserved quantity associated with translations. It is 

dL_ 



E 



n 



J] Pi 



n 



(2.57) 



which we recognise as the the total linear momentum in the direction n. Since this 
holds for all n, we conclude that ^ ■ pj is conserved. But this is very familiar. It is 
simply the conservation of total linear momentum. To summarise 



Homogeneity of Space 



Translation Invariance of L 
Conservation of Total Linear Momentum 



This statement should be intuitively clear. One point in space is much the same as any 
other. So why would a system of particles speed up to get over there, when here is just 
as good? This manifests itself as conservation of linear momentum. 

Example: Isotropy of Space 

The isotropy of space is the statement that a closed system, described by the Lagrangian 
(2.55) is invariant under rotations around an axis fi, so all are rotated by the 

same amount. To work out the corresponding conserved quantities it will suffice to 
work with the infinitesimal form of the rotations 



— Yi + ah X Yi 



(2.58) 



where a is considered infinitesimal. To see that this is indeed a rotation, you could 
calculate the length of the vector and notice that it's preserved to linear order in a. 
Then we have 



L(rj, TCi) = L(rj + an x r^, i-j + an x i-j) (2.59) 

which gives rise to the conserved quantity 
dL 

X — • (n X Ti) = X n • (r^ X Pi) = n • L (2.60) 



This is the component of the total angular momentum in the direc- 
tion fi. Since the vector n is arbitrary, wc get the result 



Figure 8: 



Isotropy of Space 



Rotational Invariance of L 

Conservation of Total Angular Momentum 
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Example: Homogeneity of Time 

What about homogeneity of time? In mathematical language, this means L is invariant 
under t — t + s or, in other words, dL/dt = 0. But we already saw earlier in this section 
that this implies H = qi{dL/dqi) — L is conserved. In the systems we're considering, 
this is simply the total energy. We see that the existence of a conserved quantity which 
we call energy can be traced to the homogeneous passage of time. Or 

Time is to Energy as Space is to Momentum 

Recall from your course on special relativity that energy and 3-momentum fit together 
to form a 4-vector which rotates under spacetime transformations. Here we see that 
the link between energy-momentum and time-space exists even in the non-relativistic 
framework of Newtonian physics. You don't have to be Einstein to see it. You just 
have to be Emmy Noether. 

Remarks: It turns out that all conservation laws in nature are related to symmetries 
through Noether's theorem. This includes the conservation of electric charge and the 
conservation of particles such as protons and neutrons (known as baryons). 

There are also discrete symmetries in Nature which don't depend on a continuous 
parameter. For example, many theories are invariant under reflection (known as parity) 
in which r.j ^ — r^. These types of symmetries do not give rise to conservation laws in 
classical physics (although they do in quantum physics). 

2.5 Applications 

Having developed all of these tools, let's now apply them to a few examples. 
2.5.1 Bead on a Rotating Hoop 

This is an example of a system with a time dependant holonomic constraint. The hoop 
is of radius a and rotates with frequency ui as shown in figure 9. The bead, of mass m, 
is threaded on the hoop and moves without friction. We want to determine its motion. 
There is a single degree of freedom ^0, the angle the bead makes with the vertical. In 
terms of Cartesian coordinates (x, y, z) the position of the bead is 

X = asinip cosLot , y = asmtlj suiLot , z = a — acosil) (2-61) 

To determine the Lagrangian in terms of the generalised coordinate ip we must substi- 
tute these expressions into the Lagrangian for the free particle. For the kinetic energy 
T we have 

T = \m{x^ + + z^) = |ma^[V^2 + u'^ sin^ V'] (2.62) 
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while the potential energy V is given by (ignoring an overall 
constant) 

V — mgz — —mgacosip (2.63) 
So, replacing x, y and z by ip, we have the Lagrangian 



L — ma' 



(2.64) 




(|)=eo t 



1 



where the effective potential is 
1 



eff 



ma" 



{—mga cos'tp 



1 2 2-2 

nTTT-a Oi) sm 



in^V) (2.65) 



Figure 9: 



We can now derive the equations of motion for the bead simply from Lagrange's equa- 
tions which read 



eflf 



(2.66) 



Let's look for stationary solutions of these equations in which the bead doesn't move 
(i.e solutions of the form ip = ip = 0). Prom the equation of motion, we must solve 
dVes/dip = to find that the bead can remain stationary at points satisfying 



g sin ip — ohP" sin cos 



(2.67) 




Figure 10: The effective potential for the bead depends on how fast the hoop is rotating 

There arc at most three such points: -f/-' = 0, i) = tt or costp = g/au^. Note that the 
first two solutions always exist, while the third stationary point is only there if the hoop 
is spinning fast enough so that o;^ > g/a. Which of these stationary points is stable 
depends on whether Ves{ip) has a local minimum (stable) or maximum (unstable). This 
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in turn depends on the value of uj. Kfr is drawn for several values of u in figure 10. 
For cu^ < (?/a, the point = at the bottom of the hoop is stable, while for uj"^ > g/a, 
the position at the bottom becomes unstable and the new solution at cos^ = gjauP' is 
the stable point. For all values of uj the bead perched at the top of the hoop ^ = tt is 
unstable. 

2.5.2 Double Pendulum 

A double pendulum is drawn in figure 11, consisting of two 
particles of mass mi and m2, connected by light rods of length 
l\ and I2. For the first particle, the kinetic energy T\ and the 
potential energy V\ are the same as for a simple pendulum 

Tx^\mxl\0\ and Vi = -mi^^i cos ^1 (2.68) 

For the second particle it's a little more involved. Consider the 
position of the second particle in the (,t, y) plane in which the 
pendulum swings (where we take the origin to be the pivot of the 
first pendulum with y increasing downwards) 

X2 = ^1 sin ^1 + ^2 sin ^2 and ^2 = cos ^1 + ^2 cos ^2 (2.69) 
Which we can substitute into the kinetic energy for the second particle 

T2 = im2(i' + Z/') 

= \m2 + ^2^2 + 2iii2 cos(^i - ^2)^1^2) (2.70) 

while the potential energy is given by 

V2 = -m^gy^ = -m^g {h cos 9i + k cos ^2) (2-71) 

The Lagrangian is given by the sum of the kinetic energies, minus the sum of the 
potential energies 

L = i(mi + m2)ll9l + \m2llel + 7^2/1^2 cos(ei - ^2)^1^2 

+(mi + m2)gli cos 9i + m25'/2 cos 62 (2.72) 

The equations of motion follow by simple calculus using Lagrange's two equations (one 
for 9i and one for 62)- The solutions to these equations are complicated. In fact, above 
a certain energy, the motion is chaotic. 
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2.5.3 Spherical Pendulum 

The spherical pendulum is allowed to rotate in three dimen- 
sions. The system has two degrees of freedom drawn in figure 
12 which cover the range 

0<e <Tr and 0<(f)<2n (2.73) 

In terms of cartesian coordinates, we have 

x — lcos(l)sm9 , y — lsijKpsmO , z — —lcos9 

Wc substitute these constraints into the Lagrangian for a free 
particle to get 

L — ^m{x^ + y'^ + i^) — mgz 

= lmf{e'^ + (f? sin^ e) + mgl cos 9 




Figure 12: 



(2.74) 



Notice that the coordinate is ignorable. From Noether's theorem, we know that the 
quantity 



J 



dL 



mP(j) sin^ 9 



(2.75) 



is constant. This is the component of angular momentum in the direction. The 
equation of motion for 9 follows from Lagrange's equations and is 



ml'^9 — mf(p sin 9 cos 9 — mgl sin 9 



(2.76) 



We can substitute cj) for the constant J in this expression to get an equation entirely in 
terms of 9 which we chose to write as 



9 



eff 



09 



where the effective potential is defined to be 



(2.77) 



(2.78) 



An important point here: we must substitute for J into the equations of motion. If 
you substitute J for directly into the Lagrangian, you will derive an equation that 
looks like the one above, but you'll get a minus sign wrong! This is because Lagrange's 
equations are derived under the assumption that 9 and are independent. 
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Figure 13: The effective potential for the spherical pendulum. 



As well as the conservation of angular momentum J, we also have dL/dt — so 
energy is conserved. This is given by 



where is a constant. In fact we can invert this equation for E to solve for 6 in terms 
of an integral 



If we succeed in writing the solution to a problem in terms of an integral like this then 
we say we've "reduced the problem to quadrature". It's kind of a cute way of saying 
we can't do the integral. But at least we have an expression for the solution that we 
can play with or, if all else fails, we can simply plot on a computer. 

Once we have an expression for 9{t) we can solve for (j){t) using the expression for 



which gives us (j) — (f){6) — (f){t). Let's get more of a handle on what these solutions 
look like. We plot the function V^fj in figure 13. For a given energy E, the particle is 
restricted to the region Ves < E (which follows from (2.79)). So from the figure we 
see that the motion is pinned between two points 6i and 62- If we draw the motion of 
the pendulum in real space, it must therefore look something like figure 14, in which 
the bob oscillates between the two extremes: 6*1 < 6* < 6*2. Note that we could make 
more progress in understanding the motion of the spherical pendulum than for the 



(2.79) 




(2.80) 
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double pendulum. The reason for this is the existence of two conservation laws for the 
spherical pendulum (energy and angular momentum) compared to just one (energy) 
for the double pendulum. 

There is a stable orbit which lies between the two 
extremal points at = 6q, corresponding to the minimum 
of VeQ. This occurs if we balance the anglular momentum 
J and the energy E just right. We can look at small oscil- 
lations around this point by expanding 9 = Oq + 59. Sub- 
stituting into the equation of motion (2.77), we have 



cflf 



(9^2 



69 + 0{59^) 



(2.81) 



0=00, 




so small oscillations about 9 = 
{d'^V^s/d9^) evaluated at ^ = 9o. 



9q have frequency o;^ 



Figure 14: 



2.5.4 Two Body Problem 

We now turn to the study of two objects interacting through a central force. The most 
famous example of this type is the gravitational interaction between two bodies in 
the solar system which leads to the eUiptic orbits of planets and the hyperbohc orbits 
of comets. Let's see how to frame this famous physics problem in the Lagrangian 
setting. We start by rewriting the Lagrangian in terms of the centre of mass R and 
the separation ri2 = ri — r2 and work with an arbitrary potential \^(|ri2|) 

L = \mj:l + \m2rl - V{\y^2\) 

= |(mi + m2)R' + |/xr?2 - ■5^(|ri2|) (2.82) 

where ^ — 17111112/ {rrii + 1112) is the reduced mass. The Lagrangian splits into a piece 
describing the centre of mass R and a piece describing the separation. This is familiar 
from Section 1.3.2. Prom now on we neglect the centre of mass piece and focus on the 
separation. We know from Noether's theorem that L = ri2 x pi2 is conserved, where 
P12 is the momentum conjugate to ri2. Since L is perpendicular to ri2. the motion of 
the orbit must lie in a plane perpendicular to L. Using polar coordinates (r, 4>) in that 
plane, the Lagrangian is 



(2.83) 



To make further progress, notice that is ignorable so, once again using Noether's 
theorem, we have the conserved quantity 



(2.84) 
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This is also conservation of angular momentum: to reduce to the Lagrangian (2.83), we 
used the fact that the direction of L is fixed; the quantity J is related to the magnitude 
of L. To figure out the motion we calculate Lagrange's equation for r from (2.83) 



Jt \'df 



dL .. -2 dV ^ 
— = fir - iir<p + — = 
or or 



(2.85) 



We can eliminate (p from this equation by writing it in terms of the constant J to get 
a differential equation for the orbit purely in terms of r, 



/^^=-|:^eff(r) 



(2.86) 



where the effective potential is given by 



Kff(r) = V{r) + 



J2 



2//r2 



(2.87) 



The last term is known as the "angular momentum barrier". Let me reiterate the 
warning of the spherical pendulum: do not substitute J = fir'^(f) directly into the 
Lagrangian - you will get a minus sign wrong! You must substitute it into the equations 
of motion. 



hyperbolic orbit 




Figure 15: The effective potential for two bodies interacting gravitationally. 



So far, you may recognise that the analysis has been rather similar to that of the 
spherical pendulum. Let's continue following that path. Since dL/dt — 0, Noether 
tells us that energy is conserved and 

£;=i//r' + Kff(r) (2.88) 
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is constant throughout the motion. We can use this fact to "reduce to quadrature", 

t-t,^^\m / (2.89) 

Up to this point the analysis is for an arbitrary potential Vir). At this point let's 
specialise to the case of two bodies interacting gravitationally with 

V{r) = ^ 2.90) 

r 

where G is Newton's constant. For this potential, the different solutions were studied 
in your Part I mechanics course where Kepler's laws were derived. The orbits fall into 
two categories: elliptic if £■ < and hyperbohc if > as shown in figure 15. 

It's worth noting the methodology we used to solve this problem. We started with 
6 degrees of freedom describing the positions of two particles. Eliminating the centre 
of mass reduced this to 3 degrees of freedom describing the separation. We then used 
conservation of the direction of L to reduce to 2 degrees of freedom (r and 0), and 
conservation of the magnitude of L to reduce to a single variable r . Finally conservation 
of E allowed us to solve the problem. You might now be getting an idea about how 
important conservation laws are to help us solve problems! 

2.5.5 Restricted Three Body Problem 

Consider three masses mi, and 777.3 interacting gravitationally. In general this prob- 
lem does not have an analytic solution and we must resort to numerical methods (i.e. 
putting it on a computer). However, suppose that 7723 <^ m\ and 7772- Then it is a good 
approximation to first solve for the motion of m\ and 7772 interacting alone, and then 
solve for the motion of 7773 in the time dependent potential set up by m\ and m2. Let's 
see how this works. 

For simplicity, let's assume ?77i and 7712 are in a circular orbit with = lot. We saw 
in the previous section that the circular orbit occurs for (9Kfj/(9r = 0, from which we 
get an expression relating the angular velocity of the orbit to the distance 



^2 ^ (2.91) 



which is a special case of Kepler's third law. Let's further assume that 7713 moves in the 
same plane as 777i and 777.2 (which is a pretty good assumption for the sun-earth-moon 
system). To solve for the motion of 777.3 in this background, we use our ability to change 
coordinates. Let's go to a frame which rotates with 777i and 7772 with the centre of 
mass at the origin. The particle 777i is a distance r^/rrii from the origin, while 7772 is a 
distance r/x/7772 from the origin. 
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Then, from the example of Section 2.2.1, the 
Lagrangian for ms in the rotating frame is 

L = |m3 [{x - uyf + (?) + uxf] - V 

where V is the gravitational potential for inter- 
acting with mi and m2 



V ^ 





y 




rn/mj 


T\ilm2 








• 




w ^ 






m2 



ri3 



r23 



(2.92) 



Figure 16: 



The separations are given by 

-"^ — [x + r n / rriif' -\- 



13 



23 



[X 



(2.93) 



Be aware that x and y are the dynamical coordinates in this system, while r is the 
fixed separation between mi and m2. The equations of motion arising from L are 



.. ^ . 2 

rrizx — zm^Luy + m^uj x — — 

ox 



mzy — —2mzOJX + mzOJ y 



dV 
dy 



(2.94) 



The full solutions to these equations are interesting and complicated. In fact, in 1889, 
Poincare studied the restricted three-body system and discovered the concept of chaos 
in dynamical systems for the first time (and, in the process, won 2,500 krona and lost 
3,500 krona). We'll be a little less ambitious here and try to find solutions of the form 
X — y — 0. This is where the third body sits stationary to the other two and the whole 
system rotates together. Physically, the centrifugal force of the third body exactly 
cancels its gravitational force. The equations we have to solve are 



rrisuj X 



msLj'^y 



dV ^ x-\-ralmi ^ x — riilmi 
— = Gmims ^ h Gm2m3 ^ 

dV 
dy 



13 



23 



(2.95) 
(2.96) 



13 



' 23 



There are five solutions to these equations. Firstly suppose that y = so that ms sits 
on the same line as mi and m2. Then we have to solve the algebraic equation 



2 ^ x + rfx/mi x-rn/m2 

UJ X — Gmi-. -. -|- Gm2i 



\x + rjji/mif 



\x 



riJ,/m2\' 



(2.97) 



In figure 17, we have plotted the LHS and RHS of this equation to demonstrate the 
three solutions, one in each of the regimes: 
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-rn/iHj 



r)i/m2 



Figure 17: The three solutions sitting on y = 0. 



r/i r/i rjj, r/i 

X < , < X < — , X > 

nil ^2 



mi mi m2 m2 

Now let's look for solutions with y ^ 0. From (2.96) we have 



Gm 



2 2 

- = cu — 



Gmi 



(2.98) 



(2.99) 



23 



13 



which we can substitute into (2.95) and, after a little algebra, we find the condition for 
solutions to be 

2 G{mi + 777,2) G(777i + 7772) 



(2.100) 



13 



23 



which means that we must have ri3 = r23 = r. There are two such points. 

In general there are five stationary points 
drawn in the figure. These are called Lagrange 
points. It turns out that Li, L2 and L3 are 
unstable, while L4 and L5 are stable as long as 
7772 is sufficiently less than mi. 

For the earth-sun system, NASA and ESA 
make use of the Lagrange points L2 and L3 to 
place satellites. There are solar observatories 
at La; sateUites such as WMAP and PLANCK 
which measure the cosmic microwave background 
radiation (the afterglow of the big bang) gather 
their data from L2. Apparently, there is a large 
collection of cosmic dust which has accumulated at L4 and L5. Other planetary sys- 
tems (e.g. the sun-jupiter and sun-mars systems) have large asteroids, known as trojans, 
trapped at their L4 and L5. 




Figure 18: The five Lagrange points. 
X marks the spots. 
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2.5.6 Purely Kinetic Lagrangians 

Often in physics, one is interested in systems with only kinetic energy and no potential 
energy. For a system with n dynamical degrees of freedom g", a = 1, . . . , n, the most 
general form of the Lagrangian with just a kinetic term is 

L^^QaMq^q' (2.101) 

The functions Qah = 9ha depend on all the generalised coordinates. Assume that 
det(g'ab) 7^ so that the inverse matrix g"-^ exists {g"'^gbc = ^"c)- ^ short exer- 

cise to show that Lagrange's equation for this system are given by 

r + C^V = (2.102) 

where 

The functions gab define a metric on the configuration space, and the equations (2.102) 
are known as the geodesic equations. They appear naturally in general relativity where 
they describe a particle moving in curved spacetime. Lagrangians of the form (2.101) 
also appear in many other areas of physics, including the condensed matter physics, the 
theory of nuclear forces and string theory. In these contexts, the systems are referred 
to as sigma models. 

2.5.7 Pcirticles in Electromagnetic Fields 

We saw from the beginning that the Lagrangian formulation works with conservative 
forces which can be written in terms of a potential. It is no good at dealing with friction 
forces which are often of the type F — —kx. But there are other velocity dependent 
forces which arise in the fundamental laws of Nature. It's a crucial fact about Nature 
that all of these can be written in Lagrangian form. Let's illustrate this in an important 
example. 

Recall that the electric field E and the magnetic field B can be written in terms of 
a vector potential A(r, t) and a scalar potential 0(r, t) 

1 (9A 

B = VxA , E = -V0--— (2.104) 

where c is the speed of hght. Let's study the Lagrangian for a particle of electric charge 
e of the form, 

L = imr^ - e (0 - ir ■ A) (2.105) 
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The momentum conjugate to r is 



— = mr+-A (2.106) 
or c 



Notice that the momentum is not simply mr; it's modified in the presence of electric 
and magnetic fields. Now we can calculate Lagrange's equations 

To disentangle this, let's work with indices a,b — 1, 2, 3 on the Cartesian coordinates, 
and rewrite the equation of motion as 

fdcP , ldAa\ , e fdAb dAa\ 



mr» = -e ^ + + _ _ (2.108) 

Now we use our definitions of the E and B fields (2.104) which, in terms of indices, 
read 

E^^-^-\^^ , B, = eJ4l (2.109) 
so the equation of motion can be written as 



mr" = eEa + -e^abBcv'' (2.110) 



or, reverting to vector notation. 



mr = e^E + -rxB^ (2-111) 

which is the Lorentz force law. 

Gauge Invariance: The scalar and vector potentials are not unique. We may make 
a change of the form 

dK 

(t>^(t^--Q^ ' A^A + cVA (2.112) 

These give the same E and B fields for any function A. This is known as a gauge 
transformation. Under this change, we have 

L^L + e^ + er- VA = L + e^ (2.113) 

ot dt 

but we know that the equations of motion remain invariant under the addition of a 
total derivative to the Lagrangian. This concept of gauge invariance underpins much 
of modern physics. 
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2.6 Small Oscillations and Stability 

"Physics is that subset of human experience which can be reduced to cou- 
pled harmonic oscillators" Michael Peskm 

Peskin doesn't say this to knock physics. He's just a fan of harmonic oscillators. And 
rightly so. By studying the simple harmonic oscillator and its relatives in ever more 
inventive ways we understand why the stars shine and why lasers shine and, thanks to 
Hawking, why even black holes shine. 

In this section we'll see one reason why the simple harmonic oscillator is so important 
to us. We will study the motion of systems close to equilibrium and see that the 
dynamics is described by n decoupled simple harmonic oscillators, each ringing at a 
different frequency. 

Let's start with a single degree of freedom x. We've already seen several examples 
where we get an equation of the form 

X = fix) (2.114) 

An equilibrium point, x — xq, of this system satisfies f{xo) — 0. This means that if we 
start with the initial conditions 

X = Xq and X = (2.115) 

then the system will stay there forever. But what if we start slightly away from x — XqI 
To analyse this, we write 

x{t) ^ xo + r]{t) (2.116) 

where r] is assumed to be small so that we can Taylor expand f{x) to find 

ri-f{xo)r) +0{rf) (2.117) 

and we neglect the terms quadratic in 7] and higher. There are two possible behaviours 
of this system 

1. f'{xo) < 0. In this case the restoring force sends us back tor) — and the solution 
is 

r]{t) = A cos{u;{t - to)) (2.118) 

where A and to are integration constants, while u'^ = —f'{xo). The system 
undergoes stable oscillations about x — xo frequency uj. 
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2. f'{xo) > 0. In this case, the force pushes us away from equihbrium and the 
solution is 



ri{t) = Ae^* + Be 



-At 



(2.119) 



where A and B are integration constants, while = f'{xo). In this case, there 
is a very special initial condition A = such that X — ^ Xq at late times. But for 
generic initial conditions, rj gets rapidly large and the approximation that t] is 
small breaks down. We say the system has a linear instability. 

Now let's generalise this discussion to n degrees of freedom with equations of motion 
of the form. 



Qi = fi{Ql,---,Qn) 



i — 1, . . . ,n 



(2.120) 



An equilibrium point g° must satisfy /i(g?, . . . , g°) = for all i = 1 . . . ,n. Consider 
small perturbations away from the equilibrium point 



Qiit) = q- + Viit) 



(2.121) 



where, again, we take the rji to be small so that we can Taylor expand the fi, and 
neglect the quadratic terms and higher. We have 



dqj 



ik=<i°k 



Vj 



(2.122) 



where the sum over j = 1, . . . , n is implicit. It's useful to write this in matrix form. 
We define the vector r] and the n x n matrix F as 



( 9h 
dqi 



dfn 
dqi 



dq„ » 



dfn , 

dqn / 



where each partial derivative in the matrix F is evaluated at 
now becomes simply 

77 = Fr) 



(2.123) 



q^. The equation 



(2.124) 



Our strategy is simple: we search for eigenvectors of F. If F were a symmetric matrix, 
it would have a complete set of orthogonal eigenvectors with real eigenvalues. Unfortu- 
nately, we can't assume that F is symmetric. Nonetheless, it is true that for equations 
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of the form (2.124) arising from physical Lagrangian systems, the eigenvalues will be 
real. We shall postpone a proof of this fact for a couple of paragraphs and continue 
under the assumption that F has real eigenvalues. In general, F will have different left 
and right eigenvectors, 

Ft^a->^lf^a , CaF^XlC a = 1, . . . , n (2.125) 

where there's no sum over a in these equations. The left and right eigenvectors satisfy 
Ca ■ Mfe = 5 ah- Note that although the eigenvectors differ, the eigenvalues for a — 
1, . . . , n are the same. Although are real for the physical systems of interest (to be 
proved shortly) they are not always positive. The most general solution to i) — Fr) is 

r,{t) = ^l^a [^ae'"' + Bae-^^'] (2.126) 

a 

where Aa and Ba are 2n integration constants. Again, we have two possibilities for 
each eigenvalue 

1. Aa < In this case ±Aa = iuja for some real number cua- The system will be stable 
in the corresponding direction t] — n^- 

2. Aa > 0. Now ±Ao are real and the system exhibits a hnear instabihty in the 
direction rj — 

The eigenvectors are called normal modes. The equilibrium point is only stable if 
Aa < for every a = 1, . . . ,n. If this is the case the system will typically oscillate 
around the equilibrium point as a linear superposition of all the normal modes, each 
at a different frequency. 

To keep things real, we can write the most general solution as 

r](t)^ J2 l^a[Aae^^' + Bae~^'^']+ J2 l^aAaCOs{u;a(t-ta)) (2.127) 

where now Aa, Ba and ta are the 2n integration constants. 
The Reality of the Eigenvalues 

Finally, let's show what we put off above: that the eigenvalues A^ are real for matrices 
F derived from a physical Lagrangian system. Consider a general Lagrangian of the 
form, 

L^lT,j{q)mj-V{q) (2.128) 
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We will require that Tij{q) is invertible for all q. Expanding about an equlibrium point 
as in (2.121), to linear order in rji the equations read 

Tijiij = -Vijrjj (2.129) 

where Tij = Tij{q^) and Vij = d^V/dqidqj, evaluated at qi = g°. Then in the ma- 
trix notation of (2.124), we have F = —T~^V. Both T^- and Vij are symmetric, but 
not necessarily simultaneously diagonalisable. This means that is not necessarily 
symmetric. Nevertheless, F does have real eigenvalues. To see this, look at 

F/li = AV =^ Vix^-X^Tn (2.130) 

So far, both /x and could be complex. We will now show that they're not. Take 
the inner product of this equation with the complex conjugate eigenvector Jx. We have 
Jx ■ Vfj, = }? Jx- Tlx. But for any symmetric matrix S, the quantity Jx ■ Sfx is real. (This 
follows from expanding /x in the complete set of real, orthogonal eigenvectors of S, each 
of which has a real eigenvalue). Therefore both Ji ■ V/x and jx ■ Tfi are both real. Since 
we have assumed that T is invertible, we know that jx ■ Tjx ^ so, from (2.130), we 
conclude that the eigenvalue is indeed real. 

2.6.1 Example: The Double Pendulum 

In section 2.5.2, we derived the Lagrangian for the double pendulum. Restricting to 
the case where the two masses are the same mi = — m and the two lengths are the 
same li — I2 — I, we derived the Lagrangian (2.72) for arbitrary oscillations 

L = ml'^9l + \ml'^6l + mt^ cos(^i - 6*2)^1^2 + '^mgl cos 6*1 + mgl cos 02 

The stable equilibrium point is clearly 61 = 62 = 0. (You could check mathematically if 
you're dubious). Let's expand for small 61 and 62- If we want to linearise the equations 
of motion for 6, then we must expand the Lagrangian to second order (so that after we 
take derivatives, there's still a 9 left standing). We have 

L ^ ml^el + \mf9l + mfOA - mgWl - \mgiel (2.131) 

where we've thrown away an irrelevant constant. From this we can use Lagrange's 
equations to derive the two linearised equations of motion 

2ml^ei + mt^h = -2mgWi 
ml% + rnl% = -mgie2 (2.132) 

Or, writing d — {9i, ^2)^, this becomes 
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Normal Mode 1 



Normal Mode 2 





Figure 19: The two normal modes of the double pendulum 




d 



d 




(2.133) 



2 

i\a 1 

We have two eigenvectors. They are 

1. {x^ = { j which has eigenvalue — —{g/l){2 — \/2). This corresponds to 
the motion shown in figure 19 for the first normal mode. 

1 



2. /U2 



-V2 



which has eigenvalue A| = —{g/l){2 + y/2). This corresponds to 



the motion shown in figure 19 for the second normal mode. 

Wc see that the frequency of the mode in which the two rods oscillate in different 
directions should be higher than that in which they oscillate together. 



2.6.2 Example: The Linear Triatomic Molecule 

Consider the molecule drawn in the figure. It's a rough 
approximation of C02- We'll only consider motion in the 
direction parallel to the molecule for each atom, in which 
case the Lagrangian for this object takes the form, 




Figure 20: 



L — ^mx\ + \Mx\ + ^mx\ — V[xi — X2) — V{x2 — x^) 



(2.134) 



The function V is some rather complicated interatomic potential. But, the point of this 
section is that if we're interested in oscillations around equilibrium, this doesn't matter. 
Assume that Xi = x° in equilibrium. By symmetry, we have — = — X3I = Tq. 
We write deviations from equilibrium as 



Xi{t) = x^, + r]i{t) 



(2.135) 
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Taylor expanding the potential about the equilibrium point, 



dV 

V{r) = F(ro) + ^ 



(r - ro) + 



r=ro 



Qj.2 



(2.136) 



r=ro 



Here the first term F (ro) is a constant and can be ignored, while the second term dV/ dr 
vanishes since we are in equilibrium. Substituting into the Lagrangian, we have 



k 



L \mf}l + \Mf]l + ^mfil - - [{jii - 7)2^ + - 773)^] 



where k — d'^V/dr^ evaluated at r = tq. Then the equations of motion are 



-k 



iV2 



V 



or, putting it in the form r) — Fr], we have 

^ —k/m k/m 



ri3 -V2 



\ 

k/M -2k /M k/M 
y k/m —k/m J 

Again, we must look for eigenvectors of F. There are three: 

Normal Mode 1 Normal Mode 2 Normal Mode 3 



(2.137) 



(2.138) 



(2.139) 



Figure 21: The three normal modes of the triatomic molecule 

1. fj, = (1, 1, 1)^ which has eigenvalue Xf — 0. But this is just an overall translation 
of the molecule. It's not an oscillation. 

2. ^2 = (IjOj"!)"^ which has eigenvalue A2 = —k/m. In this motion, the outer 
two atoms oscillate out of phase, while the middle atom remains stationary. The 
oscillation has frequency U2 = \fkfm. 

3. /Ltg = (l,-2m/M, 1)^ which has eigenvalue \% = -{k/m){l + 2m/M). This 
oscillation is a little less obvious. The two outer atoms move in the same direction, 
while the middle atom moves in the opposite direction. If M > 2m, the frequency 
of this vibration uj2, = y^— A| is less than that of the second normal mode. 
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All motions are drawn in figure 21. For small deviations from equilibrium, the most 
general motion is a superposition of all of these modes. 

r]{t) = HiiA + Bt) + n^C cos{uj2{t - t2)) + n^D cos{uj3{t - h)) (2.140) 
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